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PROCEDE DE DETERMINATION DE LA VARIATION DE LA 
PERMEABILITE RELATIVE A AU MOINS UN FLUIDE D f UN 
RESERVOIR CONTENANT DES FLUIDE S EN FONCTION DE LA 
SATURATION EN L'UN D'ENTRE EUX 

DESCRIPTION 

DOMAINE TECHNIQUE 

La presente invention est relative a 1' etude 
hydrodynamique du sous-sol en particulier dans le cadre 
de l'extraction d'hydrocarbures . Ces hydrocarbures se 
logent et se deplacent dans des formations rocheuses 
poreuses appelees reservoirs. 

De maniere inevitable, au moins un fluide 
auxiliaire se trouve present dans un reservoir 
d' hydrocarbures " en plus , de' celui que 1 ' on veut 
extraire. Le fluide a extraire peut etre du petrole ou 
du gaz, tandis que le fluide auxiliaire est 
generalernent de l'eau. Le fluide auxiliaire peut "se 
trouver naturellement en tres petite quantite dans les 
.formations rocheuses du reservoir, il s'agit d'eau 
residuelle appelee aussi eau connee, mais dans ce cas 
il est immobile. Il peut aussi .etre. .present en plus 
grande quantite • sous, forme de nappe aquifere qui se 
deplace lors de l'extraction. 

II peut enfin provenir d'un puits d' injection 
d'eau et etre utilise dans le gisement .petrolif ere pour 
refouler les hydrocarbures vers un- puits de production. 
Ce fluide auxiliaire injecte maintient.- ou remet . en 
pression 'le gisement. Cette technique d' injection de 
fluide auxiliaire est maintenant bien souvent mise en 
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place dans les premiers temps d ' exploitation du 
gisement . 

La production d'un reservoir depend des 
caracteristiques statiques du reservoir comme ses 
5 dimensions e't la nature des roches poreuses dans 
lesquelles les hydrocarbures ont pu se loger, mais 
egalement "de caracteristiques dynamiques car les 
fluides en presence se deplacent dans -le reservoir vers 
le puits lors de l 7 extraction. L'ecoulement d'un 

10 premier fluide 0 dans un reservoir nontenant un second 
fluide W est regi entre autre par la- permeabilite 
effective kefO au premier fluide O qui est le produit 
de la permeabilite absolue k du reservoir sature en 
premier fluide O par- un f acteur de correction, appele 

15 permeabilite relative krO au premier fluide - 0 en 
presence du second fluide W, cette grandeur caracterise 
la facilite avec laquelle la formation rocheuse du 
reservoir se laisse traverser par le premier fluide O 
en presence du second fluide W. 

2 0 ETAT DE LA TECHNIQUE ANTERIEURE 

De nos jours, avant de prendre une decision de 
mise en ■ " ' exploitation ' d ( ' un nouveau reservoir 
d' hydrocarbures , on 'simule- sur . ordinateur son 
comportement a 1 ' aide d'un modele dynamique de 
25 l'ecoulement des fluides dans le reservoir. Le modele 
dynamique esbt construit en trois dimensions avec une 
pluralite de 1 blocs de simulation qui sont des cubes 
dont 'la taille moyenne 'de 1' arete est de 1 ' ordre de la 
dizaine de metres. 

3 0 On a besoin de corinaitre de- la maniere la plus 

precise possible notamment pour faire la simulation la 
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variation de la permeabilite relative krO, krW a au 
moins un des fluides 0, W en fonction de la saturation 
S(W) ou S(0) en 1 ' un d'entre eux W ou 0 . La saturation 
en un fluide d'une roche est . la fraction du volume 
effectif des pores de la roche occupes par le fluide. 

Dans la pratique 1 ' un des fluides 0 est un 
hydrocarbure petrole ou gaz et 1 ' autre W est de 1 ' eau 
et c'est generalement la saturation en eau S (W) que 
1 ' on utilise comme parametre . 

Pour construire le modele dynamique de 
l'ecoulement, on . a besoin . de connaitre des 
caracteristiques statiques du reservoir et pour cela on 
. commence par construire un modele statique du reservoir 
prenant en compte des mesures geometriques du reservoir 
et -les caracteristiques geOlogiques du sous-sol, ces 
mesures peuvent etre effectuees par diagraphie. 

Le modele dynamique de l'ecoulement utilise le 
modele statique du . reservoir dans lequel on integre la 
permeabilite effective kefO de base a l'un des fluides 
0, celui que 1 ' on veut extraire, en presence d'une 
quantite residuelle. d'un autre fluide W. 

En pratique. : on utilise la permeabilite 
effective de base- au petrole. ou.au gaz en presence 

d'eau residuelle- . • : . 

Cette permeabilite effective kefO est obtenue 
de maniere classique par un essai . ,de remontee en 
pression. Dans ce type d'essai, le reservoir etant 
traverse par au moins un puits, on commence a extraire 
le fluide O avec un certain .debit, on arrete la 
production et on mesure dans, le temps la remontee en 
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pression provoquee par la fermeture du puits. On en 
deduit la permeabilite effective kefO. 

Pour faire fonctionner le modele dynamique, il 
faut 1 ' initialiser . On introduit pour cela les 
performances en debit et/ou pression. attendues pour le 
reservoir lors de son exploitation-. 

On 1' initialise egalement avec des valeurs de 
la permeabilite relative a chacun . . des fluides en 
fonction de la saturation en l'un des fluides. Ces 
valeurs de permeabilite relative doivent etre les plus 
precises possibles pour que la simulation soit fiable. 
Les valeurs de la permeabi-1 ite relative a chacun des 
fluides peuvent etre extraites de recueils de donnees 
concernant le reservoir en question, ces recueils 
existent pour la plupart les zones de prospection. Mais 
dans ce cas, il s'agit de permeabilites relatives 
moyennes . Il est meilleur au point, de vue precision 
d'effectuer des analyses de -.permeabilite, en 
laboratoire, a partir de carottes prelevees dans la 
formation du reservoir. Ces carottes sont des cylindres 
pleins, d' environ 10 centimetres de diametre, extraits 
du sous-sol. Les' mesures - de . saturation ne sont 
generaiement pas effec-tuees sur des carottes entieres 
mais sur de petits echantillons preleves dans les 
carottes. Ces carottes durant leur extraction peuvent 
subir, a cause de l'outil de coupe, des contraintes 
mecaniques irreversibles qui risquent de modifier leur 
permeabilite relative-. De plus . lors. de leur rernontee, 
elles* risquent de traverser de la boue et de s'en 
impregner, ce qui changera- egalement leur permeabilite 
relative.' Les valeurs de la permeabilite relative 
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obtenues en laboratoire different alors des valeurs de 
permeabilite relative reelles, in situ, de la formation 
du reservoir. 

De plus, les valeurs de permeabilites relatives 
determinees a partir de ces petits echantillons doivent 
etre extrapolees pour les ramener a l'echelle des blocs 
de simulation, ce qui introduit une autre source 
d ' erreurs . . 

Avec toutes ces imprecisions sur les valeurs de 
permeabilite relative, le modele dynamique risque de ne 
pas permettre de faire une representation correcte des 
ecoulernents des fluides contenus . .dans le reservoir a 
l'echelle des blocs de simulation. 

EXPOSE DE L' INVENTION 

La presente invention vise justement a 
s'affranchir des problemes mentionnes ci-dessus. La 
presente invention est un procede de determination de 
la variation de la permeabilite relative a au moins un 
fluide d'un reservoir contenant des fluides en 
mouvement en fonction- de la saturation en 1 ' un des 
fluides. Ce procede- est utilisable pour simuler, de 
maniere .plus fiable qu ' auparavant , le comportement d'un 
reservoir contenant des fluides en mouvement dont un 
des fluides est a extraire. 

Pour y parvenir le procede de determination, 
dans un reservoir contenant des fluides en mouvement, 
de la variation de la permeabilite relative a au moins 
un fluide du reservoir en fonction de la saturation du 
reservoir en 1 ' un des fluides comporte une etape de 
construction d'un modele dynamique de l'ecoulement des 
fluides dans le reservoir, destine a simuler dans le 
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temps le comportement du reservoir et d ' initialisation 
du modele dynamique avec des valeurs de base de la 
permeabilite relative a chacun des fluides en fonction 
de la saturation en l'un d'entre eux. 

Le procede est caracterise en ce qu'il comporte 
les etapes suivantes : 

a) mesure de parametres electriques du 
reservoir a 1 ' aide d'un reseau d' electrodes implante 
dans le reservoir et, a partir des mesures, obtention 
d'une reponse mesuree du reservoir a un moment donne, 
le reservoir produisant l'un des fluides, 

b) elaboration d'une reponse du reservoir de 
meme nature que la reponse mesuree, sensiblement dans 
les memes conditions, cette reponse etant simulee par 
le modele dynamique, 

c) comparaison entre la reponse mesuree et la 
reponse simulee, 

lorsque les reponses coincident, la variation 
de, la permeabilite relative audit fluide est celle que 
donne le modele dynamique dans les conditions de la 

reponse simulee, c ■ , 

en, 1' absence de coincidence, mise a jour du 
modele dynamique avec des valeurs intermediates de la 
permeabilite relative et reiteration des etapes b et c 
pour faire converger les reponses mesuree et simulee. 

La presente invention concerne egalement un 
procede de realisation d'un modele dynamique plus 
performant que ceux employes actuellement . 

Le procede de realisation d'un modele dynamique 
de l'ecoulement dans un reservoir contenant des fluides 
en mouvement, pour simuler . dans le temps le 
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comportement du reservoir. Il comporte une etape de 
construction du models dynamique de 1 ' ecoulement dans 
le reservoir et d ' initialisation avec des valeurs de 
base de la permeabilite relative a chacun des fluides 
en fonction de la saturation du reservoir en l'un des 
fluides. Il comporte egalement les etapes suivantes : 

a) mesure de parametres electriques du 
reservoir a l'aide d'un reseau d' electrodes implante 
dans le reservoir et a partir des mesures, obtention 
d'une reponse mesuree du reservoir a un moment donne, 
le reservoir produisant 1 ' un des fluides, 

b) elaboration d'une reponse du reservoir de 
meme nature que la reponse mesuree, sensiblement dans 
les memes conditions, cette reponse etant simulee par 
le modele dynamique, 

c) comparaison entre la reponse mesuree et la 
reponse simulee et tant que les reponses divergent, 
mise a jour du modele dynamique avec des valeurs 
intermediaires de la permeabilite relative et 
reiteration des etapes b et c jusqu'a ce que les 
reponses mesuree et simulee coincident. 

Il peut etre interessant, bien que les reponses 
mesuree et simulee coincident de reiterer les etapes a 
a c pour recalculer une ou plusieurs courbes de 
variation de la permeabilite relative, a differents 
moments, et pour af finer encore le f one tionnement du 
modele dynamique. 

Les valeurs de base de la permeabilite relative 
audit fluide peuvent etre obtenues a partir d' analyses 
faites sur des echantillons geologiques preleves dans 
reservoir. Selon une variante, ces valeurs peuvent etre 



SP 19524 CS 



1er depot ; 

8 

obtenues dans des recueils de donnees concernant le 
reservoir, existants auparavant. 

Les parametres electriques mesures a 1 ' aide du 
reseau d' electrodes sont des tensions ou des courants 
electriques. 

Les reponses mesuree et simulee peuvent etre 
des reponses en termes de tension ou de courant 
electriques au niveau du reseau d ' electrodes . La 
reponse simulee, en termes de tension ou de courant 
electriques, au niveau du reseau 'd' electrodes est 
obtenue a partir d'un modele de -resistivite lui-meme 
obtenu en appliquant la formule d' Archie a un- modele de 
saturation deduit du- modele d'ecoulement- dynamique. La 
reponse mesuree en termes de tension ou de courant 
electriques est obtenue directement des parametres 
electriques mesures'. . 

Dans une variante, les reponses mesuree et 
simulee peuvent etre des reponses en termes de 
resistivite apparente . La reponse mesuree en terme de 
resistivite apparente est obtenue par un calcul 
analytique simple a partir des parametres electriques 
mesures. 

La reponse simulee en termes de resistivite 
apparente est obtenue par un calcul analytique simple a 
partir d'une reponse simulee en termes de tension ou de 
courant au niveau du reseau d ' electrodes , elle meme 
obtenue a partir d'un modele de resistivite, lui-meme 
obtenu en appliquant la formule d' Archie a un modele de 
saturation deduit du modele d'ecoulement dynamique. 

Dans une autre variante, les reponses mesuree 
et simulee peuvent ' etre des reponses en termes de 
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saturation. La reponse mesuree en termes de saturation 
est un modele de saturation obtenu en appliquant la 
formule d' Archie a un modele de resistivite lui meme 
obtenu par une routine d' inversion a partir des 
parametres electriques mesures. 

La reponse simulee en termes de saturation est 
un modele de, saturation deduit du modele d'ecoulement 
dynamique . 

Le modele dynamique d'ecoulement est obtenu a 
partir d'un modele statique du reservoir contenant des 
caracteristiques statiques du reservoir. Le modele 
statique integre des donnees geologiques et 
geometriques du reservoir. Ces donnees peuvent etre 
obtenues par diagraphie 

Le modele dynamique d'ecoulement du reservoir 
integre la permeabilite effective a l'un des fluides a 
la saturation - residuelle d'un autre. La permeabilite 
effective a l'un des fluides a la saturation residuelle 
d'un autre est obtenue par un essai de remontee en 
pression effectue dans le reservoir. 

Lorsque l'un des fluides est injecte dans le 
reservoir, la reponse mesuree est obtenue a partir 
d'une mesure courante effectuee ,a ( un moment donne 
pendant 1' injection, a comparer a une mesure de base 
effectuee avant 1.' inj ect ion . Le modele dynamique 
integre pour s ' initialiser le debit et/ou la pression 
du fluide injecte. 

Lorsqu'il n'y a pas d'injection de fluide dans 
le reservoir, la reponse mesuree est obtenue a partir 
d'une mesure courante effectuee a un moment donne a 
comparer a une mesure de .base anterieure. Le modele 
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dynamique integre pour s ' initialiser la compressibilite 
du fluide autre que celui produit et le volume 
qu'occupe ce fluide. 

BREVE DESCRIPTION DES DES SINS 

La presente invention sera mieux comprise a la 
lecture de la description d'exemples de realisation 
donnes, a titre purement indicatif et nullement 
limitatif, en faisant reference aux dessins annexes sur 
lesquels : 

- la figure 1A est une coupe d'un reservoir 
auquel peuvent s'appliquer le procede; de determination 
et celui de realisation selon' 1' invention ; 

• - la figure IB montre en detail 

1 ' instrumentation d'un puits traversant le; reservoir ; 

- la figiare 2 est un organigramme exposant un 
exemple du procede de - determination et du procede de 
realisation, selon' 1' invention ; 

- les figures- 3A, 3B montrent deux variantes de 
1 ' organigramme' de la figure -2; 

' - ' - la figure ■ 4 montre plusieurs reponses en 
resistivite mesiirees du reservoir ; 

- '_'i es figures ■ 5A et 5B montrent des courbes de 
permeabilite relative de base obtenues apres avoir fait 
fonctionner le models d A ecoulement dynamique ; 

' - les figures 6A, 6B montrent les reponses en 
resistivite mesuree et simulee obtenues avec les 
courbes de permeabilite relative des figures 4A, 4B. 
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EXPOSE DETAIIiLE DE MODES DE REALISATION 
PARTICULIERS 

On va maintenant se reporter a la figure 1A qui 
montre en coupe un reservoir d'hydrocarbures auquel 
peuvent s'appliquer le procede de determination et le 
procede de realisation selon 1' invention. Le reservoir 
1 se trouve - dans des formations terrestres 2. Un puits 
3 est creuse dans les formations terrestres et traverse 
le reservoir 1. La figure IB montre schematiquement un 
detail du puits 3. Ce puits 3 peut etre un puits 
d' extraction des- hydrocarbures et/ou un puits 
d'injection et/ou un puits de mesure . On suppose que 
dans 1' exemple decrit, il s'agit d'un puits d'injection 
d'eau et qu'il est instrumented II debouche a la 
surface du sol 4.* Il peut presenter une profondeur 
variant de quelques centaines de metres a quelques 
kilometres. Le reservoir. 1 contient un premier fluide 0 
a extraire, des hydrocarbures petrole . ou gaz par 
exemple et au moins un autre fluide W qui est injecte a 
l'interieur ou qui s'y trouve naturellement . 

Le puits . 3 contient, classiquement un cuvelage 
6 et l'extremite inferieure du cuvelage est equipSe 
d'un sabot 7 de cuvelage qui le protege et facilite sa 
descente. Le puits 3 est instruments au moins avec un 
reseau 5 d' electrodes comprenant par exemple au moins 
une electrode 5-1 pour 1' injection d'un courant dans la 
formation du reservoir, au moins une electrode 5-2 pour 
le retour du courant et plusieurs electrodes 5-a, 5-b 
de mesure. Les electrodes du reseau 5 sont en contact 
electrique avec la paroi interieure du puits 3, elles 
sont generalement solidaires de la paroi exterieure du 
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cuvelage 6 mais en sont isolees electriquement dans le 
cas ou le cuvelage 6 est conducteur de 1 ' electricite . 

On peut par exemple mesurer la difference de 
potentiel V entre deux electrodes de mesure 5 -a, 5-b 
alors qu'on fait circuler dans la formation du 
reservoir un courant I connu entre 1' electrode 
d' injection 5-1 et 1' electrode de retour 5-2. La figure 
IB. represente de maniere schematique un exemple 
d'agencement des electrodes du reseau 5. Les deux 
electrodes centrales jouent le role, d' electrodes de 
mesure et les autres 5-1, 5-2 appariees par couple 
permettent 1 1 une , d'injecter un courant dans le 
reservoir et 1 ' autre d-' en • assurer le retour, pour des 
espacements differents entre ces electrodes. En fait 
les electrodes pourraient jouer al ternat ivement le role 
d'electrode d' inj ection, . d'electrode de retour ou 
d' electrode de mesure. 

Au lieu de faire des mesures de potentiel, on 
pourrait faire des mesures de courants en appliquant 
des potentiels connus aux electrodes. 

L ' instrumentation du.puits peut se faire aussi 
comme dans la .demande de,,brevet f,rangais FR-A1-2712627 
par exemple. D'autres agencements d'electrodes sont 
bien 'sur possibles. 

Pour la mise en oeuvre des - precedes selon 
1' invention, ■ .on va injecter ,1'un des fluides W (dans 
l'exemple de l'eau) dans. la. formation du reservoir 1, 
cette injection se fait dans le cuvelage 6 et le fluide 
W s'ecoule. jusqu'a la'paroi interieure du puits 3 a 
travers des orifices 6-1. -ou perforations dont est dote 
le cuvelage 6 . 
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On va decrire maintenant dans un reservoir 
contenant des fluides 0 (par exemple un hydrocarbure 
petrole ou gaz) et W (par exemple de I'eau), les etapes 
du procede de determination, de la variation de la 
permeabilite relative a au moins un fluide 0 du 
reservoir eri fonction de la saturation du reservoir en 
l'un des fluides 0 ou W. On se refere a la figure 2 qui 
represente un exemple d' organigramme du procede. Cet 
organigramme s 1 applique aussi au procede de realisation 
d'un modele de l'ecoulement des fluides dans le 
reservoir . 

On cherche a disposer d'un modele dynamique 
(bloc 20) de l'ecoulement des fluides dans le reservoir 
1." On commence par construire un modele statique du 
reservoir (bloc 10) . Ce modele statique (bloc .10) du 
reservoir prend en compte la structure du reservoir 
avec les differentes couches de formation terrestres 
composant le reservoir (bloc 12). La structure du 
reservoir est obtenue avec. des donnees geometriques et 
geologiques telles que ses dimensions, la nature des 
couches 'rocheuses qua le composent et leur porosite. 

Les donnees geometriques /et geologiques sont 
obtenues par exemple par 1 ' interpretation de mesures de 
diagraphie et/ou de donnees sismiques (bloc 11) 
effectuees dans les formations du reservoir. 

Le modele dynamique (bloc 20) du reservoir 
integre egalement la permeabilite effective kefO a l'un 
des fluides 0, celui que 1 ' on yeut extraire, a la 
saturation residuelle- en autre fluide W. La 
permeabilite effective kefO est obtenue de maniere 
classique, par un essai (bloc 13) de remontee en 
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pression dans le reservoir 1. Pour cela, le reservoir 1 
est traverse par un puits equipe anterieurement d'au 
mo ins un capteur de pression. Le capteur de pression 
porte la reference 8 sur la figure 1A, il se trouve 
dans le puits 3. Ce capteur de pression 8 est relie a 
un dispositif de pilotage et de traitement 9 des 
signaux qu'il delivre, situe par exemple en tete de 
puits . 

■ Avant de faire l'essai, le reservoir 1 produit 
1'un des fluides 0. La production est obtenue par un 
puits de production qui peut etre distinct ou non du 
puits utilise pour les mesures . Sur .la figure 1A, il 
n'y a qu'un seul puits 3 representee et on suppose qu'il 
sert a la fois pour les mesures, la production et 
1 ' inj ection . 

' L'essai consiste a arreter la production du 
reservoir en f ermant le puits de production et on 
mesure, dans le temps, la variation de la pression au 
niveau du capteur de pression 8. L ' interpretation de la 
variation de ' pression conduit a la permeabilite 
effective kefO' du reservoir au f luide 0 a la saturation 
residuelle en autre f luide W. 

On suppose egalement dans cet exemple que le 
dispositif de pilotage et de traitement 9 pilote 
egalement le reseau 5 d' electrodes et traite les 
signaux qu'il delivre. 

Ce modele dynamique (bloc 20) de l'ecoulement 
des fluides dans le reservoir doit etre initialise pour 
qu'il puisse fonctionner. Cette initialisation se fait 
en lui fournissant des valeurs de base (bloc 21) de la 
permeabilite relative (krO) b , (krW) b a chacun des 
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fluides 0, W en fonction de la saturation S(W) en l'un 
d'entre eux W. 

Ces valeurs de base de permeabilite relative 
peuvent etre extraites de la litterature (bloc 22' sur 
5 la figure 6B) . II existe des recueils de donnees 
geologiques des zones de prospection petrolif eres . Ces 
recueils peuvent prendre la forme de bases de donnees. 

On peut aussi utiliser les valeurs obtenues par 
des tests de permeabilite (bloc 22) effectues en 
10 laboratoire a partir de prelevements faits dans des 
carottes extraites du reservoir. L ' inexactitude de ces 
valeurs de base n'a pas d' importance puisque le modele 
d'ecoulement dynamique va etre mis a jour en modifiant 
ces valeurs de permeabilite relative jusqu'a ce 
15 qu'elles coincident avec les valeurs reelles. 

Pour initialiser le modele dynamique (bloc 20), 
on entre egalement les donnees (bloc 23) relatives a 
1' injection d'un des fluide, W lorsqu'elle est 
employee, c'est a dire sa pression p dans le puits 3 
20 et/ou son debit q. L' information de pression p est 
mesuree par le capteur de pression 8. 

Si 1 ' on voulait faire fonctionner le modele 
dynamique a ce stade, il faudrait egalement introduire 
le debit q' et/ou la pression p' attendues pour le 
25 fluide O a extraire dans un ou plusieurs eventuels 
puits de production. 

Mais' a ce stade, le comportement de 
l'ecoulement dans le reservoir que va simuler ce modele 
(bloc 20) va etre • eloigne de la realite a cause du 
30 manque de precision des valeurs de base de permeabilite 
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relative (krO) b , (krW) b introduites pour 

1 ' initialisation . 

Selon une caracteristique de 1' invention, on va 
recaler le modele dynamique avec des valeurs de 
permeabilite relative intermediaires, (krO)i, (krW)i de 
plus en plus approchees de celles de la realite, de 
maniere a af finer le modele dynamique de l'ecoulement, 
la simulation sera alors aussi conforme que possible 
par rapport au comportement reel du reservoir. 

Le procede selon 1 ' invention permet de 
determiner la* variation de la permeabilite relative la 
plus approchante possible de celle de la realite, cette 
derniere n'etant pas mesurable directement sur place. 

Avant de recaler le modele dynamique- (bloc 20) , 
on va effectuer une comparaison (bloc 40) en.tre une 
reponse mesuree (bloc 30) du reservoir .et une reponse 
simulee (bloc 24) du reservoir. Cette reponse simulee 
est delivree par le ; modele dynamique et elle est de 
meme nature que la reponse mesuree, sensiblement dans 
les memes conditions, c'est-a-dire qu'elle traduit 
1 • etat du reservoir a 1 ' instant de la reponse mesuree. 

La' reponse mesuree (bloc ^30) du- reservoir est 
obtenue h 'I 'aide de mesures effectuees par le reseau 5 
d' Electrodes ' implante dans le reservoir au niveau du 
puits 3 . 

Pour obtenir la reponse mesuree. (bloc 30) du 
reservoir, on commence par faire une mesure (bloc 31) 
de base dans le puits alors quMl produit l'un des 
fluides 0 et, avant 1' injection d'un autre fluide W. On 
peut par exemple mesurer la tension aux bornes des 
electrodes de mesure 5-a, 5-b alors- qu' on fait circuler 
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un courant entre les differents couples d' electrodes 5- 
1, 5-2 d' injection et de retour. On obtient une serie 
de tensions de base. Si la formation est homogene sur 
toute la profondeur du reservoir une seule mesure 
suffit en injectant du courant entre un seul couple 
d' electrodes . 

Apres la mesure de base, on injecte le fluide W 
dans le reservoir 1 via le puits 3 (bloc 32) . On 
connait le debit q d' injection et on mesure la pression 
p dans le puits comme decrit precedemment , ainsi que 
les debits q' et pression p' dans les eventuels puits 
de production. 

Pendant cette. injection, a un moment donne, on 
realise une mesure courante de tension pour un courant 
connu. Ce moment donne correspond a une certaine duree 
d'injection du fluide W (bloc 33). Avec la mesure de 
base et la mesure co.urante, on obtient une mesure utile 
(bloc 34) par sous traction . C'est cette mesure utile 
qui va donner directement ou non la reponse mesuree du 
reservoir. 

Au lieu de faire des mesures de tension, on 
pourrait faire' des mesures de courant a tension connue . 

Dans l'exemple de la figure 2, la reponse 
mesuree (bloc 30) est une reponse en saturation. 

A partir de la mesure utile (bloc 34),. par une 
routine d' inversion classique (bloc 35), on obtient un 
modele . mesure de resistivite (bloc 36) puis en 
application de la formule d' Archie (bloc 37) , un modele 
mesure de saturation (bloc 38) . La formule d'Archie, 
bien connue de 1 -homme du metier, est une formule 
empirique qui donne la resistivite d'une formation 



SP 19524 CS 



1er depot * 

18 

contenant au moins un fluide electriquement conducteur 
en fonction de la porosite de la formation, de la 
saturation en ce fluide et de la resistivite de ce 
fluide. 

C'est sur ce modele mesure de saturation que va 
s'effectuer la comparaison dans le bloc 40. 

On va voir maintenant comment obtenir la reponse 
simulee (bloc 26) du modele dynamique de l'ecoulement 
dans le reservoir. Dans 1 ' exemple de la figure 2, il 
s'agit d'une reponse simulee en' saturation. Du modele 
d'ecoulement dynamique (bloc 20) , < on - deduit directement 
un modele simule de saturation (bloc :24) . ■ 

II suffit ensuite d'effectuer une comparaison 
(bloc 40) entre la reponse mesuree (bloc 30) et la 
reponse simulee (bloc 26) sensiblement dans les memes 
conditions . 

Si les deux reponses coincident cela signifie 
que la variation de la permeabilite relative a au moins 
un des fluides O, W en fonction de la saturation S (W) 
en l'un des fluides W, contenue dans le modele 
dynamique correspond a celle- reelle du reservoir et il 
est facile' d' en avoir corinaissance (bloc 41). Le bloc 
41 delivre- des valeurs finales de permeabilite relative 
(krO)f, (krW)f.'A la premiere comparaison il s'agit des 
valeurs de- base de permeabilite relative. 

Dans cette situation, le 'modele dynamique 
obtenu peut etre considere -comme- suffisamment fiable et 
precis. Il * peut 'f onctionner : pour simuler le 
comportement du reservoir dans le futur. 

Si les deux reponses 23, 3 0 ne coincident pas, 
cela signifie qu'il faut af finer le modele dynamique 
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(bloc 20) pour que le comportement qu'il va simuler se 
rapproche encore plus de celui reel du reservoir. On va 
recaler le modele d'ecoulement dynamique avec des 
valeurs intermediaires de permeabilite relative (krO)i, 
(krW)i a. chacun des fluides O, W en fonction de la 
saturation S(W) en 1 ' un des fluides W (bloc 25). 

En faisant fonctionner le modele dynamique 
d'ecoulement avec les valeurs intermediaires de 
permeabilite relative, on elabore une nouvelle reponse 
simulee et on recommence 1 ' etape de comparaison avec la 
meme reponse mesuree (bloc 30) . Ce recalage se fait en 
changeant eventuellement plusieurs fois les valeurs de 
permeabilite . relative (bloc 25) jusqu'a ce que- les 
'reponses simulees et mesurees coincident. 

Pour ameliorer la connaissance que 1 ' on a du 
reservoir, il est preferable, lorsque la coincidence a 
ete obtenue avec une reponse mesuree correspondant un 
moment donne,. de reiterer les etapes depuis la mesure 
pour obtenir au moins une autre variation de 
permeabilite relative a , au moins un f luide dans 
d'autres conditions. . 

On effectue une .autre., mesure courante a un 
autre moment donne ; cette derniere,. conduit a une autre 
reponse mesuree, avec le modele dynamique on elabore 
une autre reponse simulee sensiblement dans les memes 
conditions que- la reponse mesuree, et 1 ' on recommence 
1' etape de comparaison jusqu'a obtenir la coincidence. 

. Cette ■ reiteration d' etapes de mesure dans le 
temps permet de rendre encore plus fiable le modele 
d'ecoulement dynamique (bloc 20). 
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Au lieu que les reponses mesuree et simulee 
soient des reponses en termes de saturation, il est 
possible que ces reponses soient des reponses en termes 
de tensions ou de courants . II est aussi possible qu'il 
n'y ait pas d' injection de l'un des fluides dans le 
reservoir. L'un des fluides W peut etre present dans le 
reservoir sous forme de nappe aquifere. On se refere a 
la figure 3A. . 

La reponse mesuree (bloc 50) est tou jours 
obtenue a 1 ' aide de mesures .effectuees par le reseau 5 
d' electrodes implante dans le puits 3. 

On commence., comme sur 1 ' organigramme de la 
figure 2, par faire une ..mesure de base dans le puits 3 
(bloc 51) . 

On e.ffectue ulterieurement une mesure. courante, 
alors que le puits 3 produit l'un des fluides O. Cette 
mesure courante . est effectuee a . un moment donne 
reference par exemple a partir de la mesure de base 
(bloc 52) . • 

. • ;On en deduit une. mesure utile en retranchant la 

mesure "de base a. la mesure courante et c'est cette 
mesure utile qui constitue la reponse mesuree (bloc 40) 
du reservoir. . •. • 

Le -modele dynamique (bloc 20) est obtenu a 
partir du modele statique (bloc 10) comme decrit a la 
figure 2. Comme il n'y a pas d' injection de fluide W 
mais que ce" dernier s'y trouve naturellement sous forme 
de nappe aquifere, on , initialise, le modele dynamique 
(bloc -20) avec la compressibili te Cw du fluide W et le 
volume Vw de la. nappe aquifere (bloc 23 ' ) . 
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II n'y a pas de changement pour son 
initialisation avec les valeurs de permeabilite 
relatives de base (kro) b/ (krw) b . 

Du modele d'ecoulement dynamique (bloc 20), on 
deduit un modele simule de saturation (bloc 24) comme 
sur la figure 2. Par application de la formule d' Archie 
(bloc 26) au modele simule de saturation, on obtient un 
modele simule de resistivite (bloc 26) du reservoir. 

A partir du modele simule de resistivite, on 
obtient • la reponse simulee du reservoir (bloc 28) en 
termes de tension ou de courant . Cela se fait par 
application d'une • formule • analytique simple, cette 
etape ne pose probleme a l'homme du metier. 

II n'y a pas de changement par rapport a la 
figure 2 pour la suite -des etapes du procede. On 
effectue la comparaison dans le bloc 40 entre la 
reponse mesuree et la J reponse simulee sensiblement dans 
les memes conditions. 

On pourrait envisager que la reponse simulee et 
la reponse mesuree soient ■ des reponses en termes de 
resistivite apparente : C'est cette variante qu'illustre 
la figure 3B. Cette figure est- comparable d'une part a 
la figure 2 et d' autre part a la figure 3A. Les mesures 
s'effectuent comme sur la figure 2 avec injection de 

fluide W. * 

Avec la mesure de base (bloc 31) et la mesure 
courante (bloc 33) on obtient la mesure utile (bloc 
34) En appliquaht une formule analytique simple a la 
mesure utile,- on - obtient la reponse mesuree du 
reservoir en termes de resistivite apparente (bloc 60). 
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L'obtention de la reponse simulee du modele 
dynamique en termes de resistivite apparente (bloc61) 
est issue de la reponse simulee du reseau d' electrodes 
(bloc 28 sur la figure 3A) a laquelle on applique une 
formule analytique simple. 

Une autre difference dans l'organigramme de la 
figure 3B par rapport a ceux des figures 2 ou 3 A est 
que les valeurs de base de permeabilite relative 
(krO) b , (krW) b sont maintenant issues de la litterature 
au lieu d'essais sur des echantillons (bloc 22'). 

La figure 4 montre cinq reponses mesurees en 
termes de resistivite apparente du reservoir a 
differents moments distincts,. apres les durees 
d' injection de fluide W suivantes : 

Tl = 4 heures 

T2 ■ = 8 heures 

T3 = 12 heures 

T4 = 16 heures 

T5 = 2 0 heures. 

On suppose que ces reponses ont ete obtenues 
dans un reservoir d' environ 15 metres d'epaisseur, 
equipe d'un .reseau de 14 electrodes regulierement 
reparties sur toute sa. profondeur. Ce reservoir a une 
resistivite de base d'environ 500 ohm.m. On injecte 
dans le reservoir un fluide W avec un debit d' injection 
d'environ 9, 2. 10" 5 m 3 /s. En- faisant fonctionner le 
modele dynamique, on s' attend a ce .que les reponses 
simulees en termes de resistivite apparente qu'il 
produit dans des conditions corr espondant aux. 
differents moments coincident avec les reponses 
mesurees . 
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Les figures 5A et 5B representent 
respect ivement un jeu de courbes de base de la 
permeabilite relative krW en fluide W (identifiee avec 
des A) et de la permeabilite relative krO en fluide 0 
(identifiee avec des •) en fonction de la saturation 
S(W) en fluide W dans un reservoir. Les courbes en 
traits pleins correspondent aux valeurs de 
permeabilites relatives de base. Ce sont ces valeurs 
qui servent a initialiser le modele d'ecoulement 
dynamique. Les courbes en pointilles correspondent aux 
valeurs intermediates introduites lors du recalage, 
elles bnt conduit a la coincidence entre les reponses 
mesuree et simulee du reservoir pour une mesure a un 
instant donne. 

Sur les figures 6A, 6B a examiner en parallele 
avec les figures 5A, 5B, on a* represents des reponses 
simulees en termes de resistivite apparente du 
reservoir en fonction de l'ecart E entre les electrodes 
contribuant a la circulation du courant. La reponse 
simulee en traits pleins est la premiere reponse 
simulee obtenue, • en mettant en oeuvre le procede de 
determination selon 1' invention, a partir des valeurs 
de base de permeabili te relative. La reponse simulee en 
traits pointilles est la reponse simulee qui coincide 
avec la reponse mesuree, elle est obtenue apres 
plusieurs mises a jour du modele d'ecoulement 
dynamique. Apres la modification des valeurs de la 
permeabilite relative krW au fluide W de la figure 5A, 
la premiere reponse -simulee de la figure 6A a re joint 
la reponse mesuree. De la meme maniere, apres la 
modification des valeurs de la permeabilite relative 
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krO au fluide 0 de la figure 5B, la premiere reponse 
simulee de la figure 6B a re joint la reponse mesuree. 

Bien que certains modes de realisation de la 
presente invention aient ete representes et decrits de 
fagon detaillee, on comprendra que differents 
changements et' modifications puissent etre apportes 
sans sortir du cadre de 1' invention. 
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REVENDICATIONS 

1. Procede de determination, dans un reservoir 
(1) contenant des fluides (W, 0) en mouvement, de la 
5 variation de la permeabilite relative (krO, krW) a au 
moins un fluide du reservoir (1) en fonction de la 
saturation du reservoir (1) en 1 ' un des fluides (W) , 
comportant une etape de construction d'un modele 
dynamique (20) de l'ecoulement des fluides dans le 

10 reservoir (1), destine a simuler dans le temps le 
comportement du reservoir et d ' initialisation avec des 
valeurs de base (krO) b / (krW) b de la permeabilite 
relative a chacun des fluides (W,0) en fonction de la 
saturation en 1'un d'entre eux (W) , caracterise en ce 

15 qu'il comporte les etapes suivantes : 

a) mesure de parametres electriques du 
reservoir (1) a 1 ' aide d'un reseau d'electrodes (5) 
implante dans le reservoir (1) et a partir des mesures 
obtention d'une reponse mesuree (3 0) du reservoir a un 

2 0 moment donne, le reservoir produisant l'un des fluides 
(O) , 

b) elaboration d'une reponse (24) du reservoir 
de meme nature que la reponse mesuree, sensiblement 
dans les memes conditions, cette reponse etant simulee 

2 5 par le modele dynamique (20) , 

c) comparaison (40) entre la reponse mesuree 
(30) et la reponse simulee (24), 

lorsque les reponses coincident, la variation 
de la permeabilite relative audit fluide est celle que 
30 donne le modele dynamique (20) dans les conditions de 
la reponse simulee (24), 
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en 1' absence de coincidence, mise a jour du 
modele dynamique (20) avec des valeurs intermediaires 
(krO)i, (krW)i de la permeabilite relative et 
reiteration des etapes b et c pour faire converger les 
5 reponses mesuree et simulee. 

2. Procede de realisation d'un modele dynamique 
de l'ecoulement dans un reservoir (1) contenant des 
fluides (0,W) en mouvement, pour simuler dans le temps 

10 le comportement du reservoir.,, comport ant une etape de 
construction du modele. dynamique (2 0) de l'ecoulement 
des .fluides dans *le reservoir (1) et d ' initialisation 
avec des valeurs ; .de. , base (krO) b '/ (krW) b de la 
permeabilite relative a chacun des fluides (0,W) en 

15 fonction de la- saturation -du reservoir en 1 ' un des 
fluides (W) , caracterise en ce qu'il comporte les 
etapes suivant.es :-<...',* 

a) mesure • de parametres electriques du 
reservoir (1) a 1 ' aide d'un reseau (5) d' electrodes 

2 0 implante dans le reservoir et obtention d'une reponse 
mesiliree (30). du • reservoir . a. un moment donne, le 
'reservoir produisant. l'un .des fluides, • 

b) ; elaboration d'une reponse (24) du reservoir 
de meme nature que la reponse mesuree,, sensiblement 

25 dans les merries conditions, cette reponse etant simulee 
p'ar, le modele dynamique, 

c) comparaison (40) entre la reponse mesuree et 
•la reponse • simulee, et tant que les reponses divergent, 
mise <a. jour du . modele dynamique (20) avec des valeurs 

30 intermediaires de la permeabilite - (krO)i, (krW)i et 
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reiteration des etapes b et c jusqu'a ce que les 
reponses mesuree et simulee coincident. 

-3. Procede selon l'une des revendications 1 ou 

2, dans lequel, en cas de coincidence, on reitere au 
moins une fois les etapes a a c, la reponse mesuree 

(30) etant obtenue a un autre moment donne. 

4. Procede selon l'une des revendications 1 a 

3, caracterise en ce que les valeurs de base (krO) b , 
(krW) b de la • permeabili te relative audit fluide sont 
obtenues a partir (22) d'analyses faites sur des 
echantillons geologiques preleves dans reservoir. 

5. Procede selon l'une des revendications 1 a 
3, caracterise en ce que- les valeurs de base (krO) b / 
(krW) b de la permeabilite relative sont obtenues dans 
des recueils de donnees (22') concernant le reservoir, 
existants auparavant . 

6. Procede selon l'une des revendications 1 a 
5, caracterise 1 en ce que les pararnetres electriques 
mesures a 1'aide du reseau (5)- d'eTectrodes sont des 
tensions ou- des courants electriques'. 

7 . Procede selon l'une des revendications 1 a 6, 
caracterise en ce que les reponses mesuree (50) et 
simulee (28) sont des reponses en termes de tension ou 
de courant electriques au niveau du reseau 
d ' electrodes . ' -' • 
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8. Procede selon l'une des revendications 1 a 
6, caracterise en ce que les reponses mesuree (60) et 
simulee (61) sont des reponses en termes de resistivite 
electrique apparente . 

9. Procede selon l'une des revendications 1 a 
6, caracterise en ce que les reponses mesuree (20) et 
simulee (24) sont des reponses en termes de saturation. 

.10. Procede selon. la - revendication 7, 
caracterise en ce que la reponse simulee (28) en terme 
de tension ou de courant electriques.au niveau du 
reseau " (5) d'eleetrodes est obtenue a partir d'un 
modele . de. resistivite ■ ; (27), lui-meme obtenu en 
appliquant la formule d' Archie (26) a un modele de 
saturation (24) deduit du modele d'ecoulement dynamique 
(20) . . 

. • 11. Procede selon.- la revendication 7 , 
caracterise en ce '.que la reponse mesuree (50) en termes 
de tension de courant electriques est obtenue 
directement des parametres electriques mesures. 

•.12-i,'. ■ Procede ■■ selon la , revendication 8, 
caracterise en ce que la reponse mesuree (60) en termes 
de resistivite apparente- est obtenue par un calcul 
analytique simple a partir ..des parametres electriques 
mesures . 

'13. Procede selon la revendications 8, 
caracterise- en ce que la . reponse simulee (61) en termes 
de resistivite apparente est obtenue par un calcul 
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analytique simple a partir d'une reponse simulee en 
termes de tension ou de courant electriques au niveau 
du reseau (5) d ' electrodes , elle-meme obtenue a partir 
d'un modele de resistivite (27) , lui-meme obtenu en 
appliquant la formule d' Archie (26) a un modele de 
saturation (24) deduit du modele d'ecoulement dynamique 
(20). 

14. Procede selon la revendication 9, 
caracterise en ce que la reponse mesuree (30) en termes 
de saturation est un modele de saturation obtenu en 
appliquant la formule d' Archie (37) a un modele de 
resistivite (36) lui' meme obtenu par une routine 
d' inversion a partir des parametres electriques 
mesures . 

15. Procede selon la revendication 9, 
caracterise en ce que la reponse simulee (24) en termes 
de saturation est un modele. de saturation deduit du 
modele d'ecoulement dynamique (20). 

16. Procede selon l'une des revendications 1 a 
15, caracterise en ce que le modele dynamique (20) 
d'ecoulement est ■ obtenu a partir d'un- modele statique 
(10) du reservoir- contenant des caracteris t iques 
statiques du reservoir. 

17. Procede selon la revendication 16, 
caracterise en ce que le modele statique (10) integre 
des donnees (11, 12) geologiques et geometriques du 
reservoir ( 1 ) . - - 
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18. Procede selon l'une des revendications 1 a 
17, caracterise en ce que le modele 'dynamique (20) 
d'ecoulement ' du reservoir integre la permeabilite 
effective (Kefo) a 1'un des fluides (O) a la saturation 
residuelle d'un autre (W) . 

19. Procede selon la revendication 18, 
caracterise en ce que la permeabilite effective (Kefo) 
a l'un des fluides (0) a la saturation residuelle d'un 
autre (W) est obtenue par un essai de remontee en 
pression (13) effectue dans le reservoir. 

20. Procede selon l'une des revendications 1 a 
19, lorsque l'un (W) des fluides est injecte dans le 
reservoir (1), la reponse mesuree (30) est obtenue a 
partir d'une mesure courante (33) effectuee a un moment 
donne pendant 1' injection a comparer a une mesure de 
base (31) effectuee avant 1' injection. 

21. Procede selon la revendication 20, 
caracterise en ce que le modele dynamique integre pour 
s ' initialiser le debit (q) et/ou la pression (p) du 
fluide (W) injecte. 

22. Procede selon l'une des revendications 1 a 
19, caracterise en ce qu'en absence d' injection de 
fluide dans le reservoir (1), la reponse mesuree (50) 
est obtenue a partir d'une mesure courante (52) 
effectuee a un moment donne, a comparer a une mesure de 
base (51) anterieure. 
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23. Procede selon la revendication 22, 
caracterise en ce que le modele dynamique (20) integre 
pour s ' initialiser la compressibilite (C w ) du fluide 
5 (W) autre que celui produit et le volume (V w ) qu'occupe 
ce fluide (W) . 



SP 19524 CS 



1 er depot "? 



1/6 




1er depot 



2/6 



PUITS 



13 

S 

Essai de 
remontee 
en pression 



Mesure de base 



Injection fluide 



1 



^31 

^32 
.33 



Mesure courante < 



34 



Mesure utile 



Routine 
d'inversion 



.35 



36 



Modele de resistivite 



I 



Archie 



.37 



Modele de saturation 

7 



30 



22 



Donnees sismiques 



i 



Couches 



kofO 



.12 



Modele statique 



21 



Carotte 



(Kro)b 
(Krw) b 



(Kro)i 
(Krw)i 



1 4 



Modele 
d'ecoulement 
^ dynamique 



XZ-25 




20 



24 



(Kro) { 
(Krw) f 
I 

I I 



Modele de 
saturation simule 



Fig 0 2 



1er depot 



3/6 



PUITS 



13 



Essai de 
remontee 
en pression 



Mesure de base 



.51 



52 



22 



21 



Carotte 



(Kro) b 
(Krw) b 



Mesure courante «S 



Mesure utile 



X^ 



50 



(Kro)i 
(Krw)i 



25 



Donnees sismiques 



11 



Couches 



.12 



Modele statique 



Y Y 



.10 
23 



Mpdele 

d'ecoulement 
dynamique 



33 





Cw 




et 




Vw 



Modele de 
saturation simul€ 



20 

^-24 



Archie 




.25 



Modele de resistivite 



OUI 
1 

(Kro) f . 
(Krw) f 



.41 



Reponse simulee 
reseau d'electrodes 



Fig 0 3A 



.27 



28 



L _-_ J 



1er depot 



4/6 



3 



purrs 



13 

s 

Essai de 
remontee 
en pression 



Mesure de base 



,31 



32 



Injection fluide 



22 

A. 



Donriees sismiques 



.11 



Couches 



-12 



Modele statique 



21 



Iitterature 



(Kro) b 
(Krw)b 



Mesure courante «S — | 



v r- 



34 



Mesure utile 



60 



(Kro)i 
(Krw)i 



25 Z7 



i t 



.10 
23 



Modele 

d'ecoulement 
dynamique 



<- 


p ♦/ 

et/ou 




q 



20 



Modele de saturation 



,24 



Archie 



_26 
27 



Modele de resistivite 



Resistivites apparentes 




Reponse simulee 
reseau d' electrodes 



28 



61 



OUI 

, y , 

(Kro) f 
(Knr), 



Resistivites apparentes 
simulees 



Figo 3B 



.41 



I 



r 



1er depot 



6/6 





IWl 



J IKSTITUT 
NATIONAL DE 
LA PROPRUTE 
IMDU8TRIELIE 



regue le 26/11/01 I 

iC^EWET ©'OfiWEMTflOlM 
CEOTFOCATT O'yTOLDTE 

Code de la propriete intellectuelle - Livre VI 



N° 11235*02 



DEPARTEMENT DES BREVETS 

26 bis. rue de Saint Petersbourg 
75800 Paris Cedex 08 

Telephone : 01 53 04 53 04 Telecopie : 01 42 93 59 30 



DESIGNATION D'INVENTEUR(S) Paga N° .1. . / 1. . 

(Si le demandeur n'est pas I'inventeur ou I'unique inventeur) 



Cet imprime est a remplir lisiblement a l f encre noire 



D3 113 V.' /260899 



V s references pour ce dossier 

(facultatij) 



SP 19524/CS 21.0910 



W° D'EWREGISTREftflENT NATIONAL 



TITRE DE L'INVENTION (200 caracteres ou espaces maximum) 

PROCEDE DE DETERMINATION DE LA VARIATION DE LA PERMEABILITE RELATIVE A AU 
MOINS UN FLUIDE D'UN RESERVOIR CONTENANT DES FLUIDES EN FONCTION DE LA 
SATURATION EN L'UN D'ENTRE EUX. 



LE(S) DEiVJAMDEUR(S) : 

SERVICES PETROLIERS SCHLUMBERGER 
42 rue St Dominique 
75007 PARIS 





(DESIGN E(NT) EN TAWT QU'5WVENTEUR(S) : (Indlquez en haut a drolte « Page N° 1/1 » S'il y a plus de trois inventeurs, 
utilisez un formulaire identique et numerotez chaque page en indiquant le nombre total de pages). 


Nom 


DELHOMME 


Prenoms 


Jean-Pierre 


Adresse 


Rue 


20 avenue des Mimosas 


Code postal et ville 


92 1 00 BOULOGNE BILLANCOURT 


Societe d'appartenance (facullatif) 




Nom 


MANIN j 


Prenoms 


Yves 


Adresse 


Rue 


4, Grand Place 


Code postal et ville 


92350 LE PLESSIS ROBINSON. 


Societe d'appartenance (facullatif) 




Nom 




Prenoms 




Adresse 


Rue 




Code postal et ville 


1 


Societe d'appartenance (facultatij) 




DATE ET SIGNATURE(S) 
DU (DES) DERflAWDEUR(S) 
OU DU SV3ANDATAIRE 
(Nom et qualite du signataire) 

PARIS LE 8 NOVEMBRE 2001 
G.POUL1N CPI 990200 





La loi n°78-17 du 6 janvier 1978 relative a I'informatique, aux fichiers et aux libertes s'applique aux reponses faites a ce formulaire. 
Elle garantit un droit d'acces et de rectification pour les donnees vous concernant aupres de I'lNPI. 



